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Abstract 

Effective  properties  of  plain  stratified  bianisotropic  structures  are  found  on  the  basis  of  4x4  matrix 
method  without  restriction  on  the  ratio  of  the  period  of  structure  and  the  wavelength  of  incident 
radiation.  The  effective  permittivity,  magnetic  permeability  and  magnetoelectric  dyadics  are 
obtained  for  isotropic  two-layer  and  anisotropic  helicoidal  structures  in  resonant  case. 

1.  Introduction 

The  problem  of  effective  properties  of  stratified  media  is  known  in  optics  for  a long  time  [1].  The 
effective  properties  of  an  inhomogeneous  plane  stratified  structure  are  defined  as  the  properties  of  a 
certain  homogeneous  layer  which  cannot  be  distinguished  from  the  investigated  inhomogeneous 
structure  by  optical  means.  As  a rule,  periodic  structures  with  large  number  of  periods  are  considered 
in  long  wave  approximation,  i.e.  in  the  case  when  the  length  of  electromagnetic  wave  is  much  greater 
than  the  period  of  the  structure.  Such  structures  exhibit  anisotropic,  chiral  and,  generally,  bianisotropic 
properties  [1-3]. 

In  case  when  the  wavelength  of  light  is  approximately  equal  to  the  period  of  structure,  the 
resonant  interaction  of  light  with  the  structure  is  observed.  This  problem  is  solved  by  various  methods, 
for  example,  by  the  method  of  coupling  waves  or  Bloch  functions  formalism  [4].  The  method  of 
effective  properties  was  not  earlier  used  in  this  case,  being  limited  in  application  by  media  with  short 
period. 

The  purpose  of  the  present  work  is  to  find  the  effective  properties  of  periodic  bianisotropic 
structures  having  the  period  approximately  equal  to  the  wavelength  of  propagating  light,  and  also  non- 
periodic structures  made  of  a few  layers.  For  solution  of  the  given  problem  we  apply  a modification  of 
the  Berreman  4x4  matrix  method. 

2.  Method  of  4x4  Matrices;  Matrix  of  Material  Properties 

Consider  a medium  consisting  of  bianisotropic  layers  parallel  to  XV  plane.  A plane  electromagnetic 
monochromatic  wave  propagates  in  the  medium.  Its  wave  vector  k is  parallel  to  XZ  plane.  Maxwell 
equation  can  be  transformed  to  the  following  4x4  differential  matrix  equation: 

^ = -ikoG^e,fi,dJ'jg,.  (1) 

where  kg  is  the  wave  number  in  vacuum;  the  four  dyadics  of  permittivity  e , permeability  fi  and 

magnetoelectric  d , P describe  the  relations  between  electric  and  magnetic  fields;  a four-component 

vector  g = {E^,-Ey,  is  composed  of  tangential  components  of  the  fields.  Here  the  4x4  matrix 

of  material  properties  (MMP)  G is  determined  by  local  properties  of  the  medium  and  the  incidence 
angle.  It  is  constructed  on  the  basis  of  the  four  constitutive  dyadics  and  allows  in  general  way  taking 
into  account  bianisotropic  properties  of  the  medium.  MMP  consists  of  the  sum  of  three  matrices  which 
are  proportional  to  zero,  first  and  second  degrees  of  the  tangential  component  of  the  wave  vector  : 

^ ^ ^ kl  ^ 

G-Gq  +-^Gi  + — G2. 

^0  Kq 


(2) 
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For  a homogeneous  medium  MMP  does  not  depend  on  z,  and  the  solution  of  the  matrix 
equation  (1)  in  this  case  is  a superposition  of  eigen  waves: 

8 = = (3) 

where  Cj  is  the  amplitudes  of  the  waves  corresponding  to  the  eigen  vectors  of  matrix  G.  The  eigen 
numbers  of  this  matrix  = k^/k^  are  the  solutions  of  the  dispersion  equation  det(G  - n^I) , where  I 

is  the  unit  matrix.  As  follows  from  (3),  there  are  four  eigen  waves  with  various  polarizations 
determined  by  the  fourth  order  equation  with  respect  to  n^- . 

If  the  structure  contains  several  homogeneous  layers,  the  investigation  of  propagation  of 
electromagnetic  waves  turns  into  the  solution  of  a boundary  problem.  The  continuity  conditions  for 
tangential  electrical  and  magnetic  field  components  must  be  satisfied  on  the  boundaries  of  adjoining 
layers,  that  is  identical  to  the  continuity  of  the  four-component  vector  g.  Let  a wave  described  by 
vector  g(z„)  is  incident  on  m-th  layer.  The  vector  g(z„+,)  of  the  wave  transmitted  through  the  layer 
is  determined  by  the  matrix  relation: 

g(z„,,)=S<">g(z„),  (4) 

where  the  transformation  matrix  for  m-th  layer  of  thickness  Az„  has  the  form: 

S('”)=exp[-//:oG(z,)Az^].  (5) 

Here  the  exponential  function  is  applied  to  a 4x4matrix.  It  is  known  from  the  theory  of  matrices  that 

for  arbitrary  matrix  A and  function/the  following  relation  is  satisfied:  ' 

/(QAQ)  = Qf(A)Q,  (6) 

where  Q is  the  matrix  inverse  for  Q.  Let's  take  Q as  a matrix  which  are  formed  by  columns  of  eigen  ’ 

vectors  of  G:  Qjj^  =(q/^)j . Then  the  lines  of  Q are  the  vectors  complimentary  to  q^ 

(Qj’Qjk  )■  We  substitute  MMP  in  (5)  as:  G = QQGQQ.  Taking  into  account  that 
q;Gq,=q  y”z*q*  transformation  matrix  of  m-th  layer  can  by  presented  as: 

. (7) 

u 

For  a system  consisting  of  p-l  layers  (a  wave  is  incident  from  medium  n = 0 , passes  p boundaries 
and  propagates  to  the  medium  m = p)  the  resulting  transformation  matrix  is  the  product  of 
transformation  matrices  of  separate  layers: 

S = (8) 

For  inhomogeneous  medium  MMP  is  a continuous  function  of  z,  and  the  resulting  transformation 

matrix  can  be  found  as  the  product  of  transformation  matrices  of  indefinitesimal  layers,  into  which  the 
medium  is  divided.  When  the  thickness  of  these  layers  tends  to  zero  we  obtain  the  transformation 
matrix  as  a multiplicative  integral: 

S = lim  n = f g-''*oG(z)rfz  ^ (9) 

J 

where  multiplicative  integral  is  designated  as: 

[A(z)*  = lim  1~T  [l -I- A(z)Az]=  lim  nA(z  )^  . 

J J 

The  same  result  can  be  obtained  by  separating  the  variables  and  integrating  the  basic  equation  (1). 
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3.  Effective  Matrix  of  Material  Properties 


Arbitrary  stratified  bianisotropic  structure  accordingly  to  (8)  can  be  described  by  a transformation 
matrix  S,  which  contains  the  complete  information  on  all  optical  characteristics  of  the  structure  as 
whole.  This  allows  replacing  such  structure  by  a certain  homogeneous  layer  with  effective  properties 
so  that  the  transformation  matrix  does  not  changes.  Thus,  in  order  to  find  the  effective  properties  of 
the  structure  it  is  necessary  at  first  to  find  the  transformation  matrix  from  the  properties  of  the 
inhomogeneous  structure,  then  to  solve  the  inverse  problem  and  find  effective  MMP  assuming  the 
structure  to  be  a homogeneous  layer.  From  MMP  it  is  possible  to  obtain  permittivity,  permeability  and 
magnetoelectric  dyadics. 

Equating  the  transmission  matrices  of  a homogeneous  layer  (5)  and  an  inhomogeneous  structure 
with  arbitrary  function  G(z)  (9)  and  taking  the  logarithm  we  obtain  the  effective  MMP: 


UqL 


-In 


G(z)dz 


(10) 


where  L is  the  whole  thickness  of  the  structure.  The  obtained  expression  is  exact.  It  allows  calculating 
effective  properties  of  structures  both  with  large  number  of  layers  and  few  layers,  for  example,  two  or 
three-layer  structures.  For  periodic  medium  it  is  enough  to  take  integral  over  the  period  of  the 
structure. 

It  is  known  that  the  logarithmic  function  is  multivalued,  in  particular,  for  scalars  it  is  determined 
up  to  litim , where  m is  an  integer  number.  The  multyvaluedness  of  logarithmic  function  is  the  result 
of  the  periodicity  of  exponential  function.  For  a matrix  of  dimension  sxs  there  exist  ^ independent 
integer  numbers  m^,K  , varying  which  it  is  possible  to  obtain  new  matrices  so  that  the  exponential 
function  upon  these  matrices  is  a constant: 


A'  = A-l-2raQMQ, 


M = 


WJ] 

0 

0 


0 0 ' 

O 0 
0 


(11) 


where  the  matrix  Q is  made  of  the  columns  of  eigen  vectors  of  the  matrix  A.  Thus,  the  logarithm  of  a 
matrix  is  a set  of  counted  number  of  matrices.  Ambiguity  of  determination  of  MMP  follows  from  this 
fact.  Determination  of  MMP  is  ambiguous  not  only  when  we  find  effective  MMP  for  inhomogeneous 
structure  but  also  for  usual  homogeneous  medium.  Really,  let  we  transform  G — > G'  so  that 
exp(-ikoG)  = exp(-i'A:oG')  ■ Then  the  matrix  S (5)  and,  therefore,  and  reflection  and  transmission 
coefficients  will  not  change.  Hence,  the  determination  of  constitutive  dyadics  of  medium  is  ambiguous 
at  fixed  optical  characteristics. 

The  effective  MMP  for  an  inhomogeneous  medium  is  a complex  function  of  , and,  generally, 
it  cannot  be  presented  in  the  form  of  (2)  when  G,-  is  independent  on  k^.  Therefore  the  effective 

dyadics  e,  p.,  d and  p at  any  splitting  of  G into  three  summands  depend  on  k^  and  on  the  direction 
of  wave  propagation.  Hence,  such  medium  exhibit  spatial  dispersion. 


4.  Examples  of  Stratified  Structures 

Consider  two  most  typical  examples  demonstrating  application  of  the  general  approach  of  calculation 
of  effective  properties  of  layered  bianisotropic  media. 


4.1  Two-layer  structure 

We  consider  a structure  which  consists  of  two  dielectric  layers  with  permittivities  e,2=e±Ae/2 
( Ae  = e ) and  thickness  d^  2 ■ Non-zero  components  of  the  effective  dyadics  of  such  structure  can  be 
written  as: 
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The  propagation  constants  of  p and  s polarized  waves  in  such  medium  =±C7-s/£  (l 

describe  attenuating  waves.  The  positive  direction  wave  propagating  in  such  medium  attenuates  and 
transfers  its  energy  to  the  negative  direction  wave,  i.e.  Bragg  reflection  of  light  from  periodic  structure 
takes  place.  The  penetration  depth  of  the  wave  in  the  structure  is  determined  by  the  parameter  v and 

equals  nj ■ The  effective  medium  exhibits  the  property  of  strong  spatial  dispersion, 

since  constitutive  dyadics  are  functions  of  k^,  i.e.  functions  of  the  wave  vector.  This  is  consisted  with 
the  fact  that  the  contribution  of  spatial  dispersion  is  essential  when  the  relation  between  the  period  of 
stracture  and  the  wavelength  is  close  to  unit. 


4.2  Helicoidal  structure 


A Helicoidal  structure  represents  a twisted  medium  in  which  there  is  a rotation  of  anisotropy  axes 
around  Z axis  at  displacement  along  its  direction.  The  effective  properties  of  such  medium  correspond 
to  a chiral  medium  with 

««  = =-Pyy-  'V A)  - (13) 

where  L is  the  period  of  the  structure;  X is  the  wavelength  of  the  incident  radiation. 

5.  Conclusion 

On  the  basis  of  4x4  matrix  method  the  exact  equation  for  effective  MMP  of  arbitrary  bianisotropic 
structure  is  obtained  without  restriction  on  the  relation  between  the  wavelength  and  the  thicknesses  of 
the  layers  constituting  the  structure. 

The  effective  constitutive  dyadics  for  a two-layer  isotropic  dielectric  structure  are  found  in  the 
resonant  case.  It  is  shown  that  the  effective  magnetoelectric  dyadics  contain  off-diagonal  components 
which  were  absent  for  the  isotropic  layers  of  the  structure.  At  this  conditions  the  permittivity  and 
permeability  become  complex  quantities,  that  results  in  strong  Bragg  reflection  of  a wave  from  such 
stmcture.  The  effective  constitutive  dyadics  for  helicoidal  medium  are  found.  It  is  shown  that  such 
medium  is  chiral. 

MMP  and  constitutive  dyadics  of  any  medium  are  determined  ambiguously.  The  transformation 
for  MMP  which  completely  conserves  the  optical  characteristic  of  the  layered  structure  is  found. 
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